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Abstract 

The recent measurement on the decay constant of Dg shows a discrepancy between theory and 
experiment. We study the leptonic and semileptonic decays of D and Dg simultaneously within the 
standard model by employing a lightfront quark model. There is space by tuning phenomenological 
parameters which can explain the "/d^ puzzle" and do not contradict other experiments on the 
semileptonic decays. We also investigate the leptonic decays of D and Dg with a new physics 
scenario, unparticle physics. The unparticle effects induce a constructive interference with the 
standard model contribution. The nontrivial phase in unparticle physics could produce direct CP 
violation which may distinguish it from other new physics scenarios. 
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I. INTRODUCTION 



Thanks to the recent experimental improvements from B factories, BES, CLEO-c and 
hadron colUde., ch^„> ph^^cs ente.= ,nto a "golden t.„>e" Q. M.ny new charmed reso- 
nances are observed, such as Dgj and X, Y, Z states. They open a new window to study 
nonperturbative QCD. The last neutral meson mixing, D'^ — mixing is observed at about 
1% level. This is a typical flavor-changing-neutral-current (FCNC) transition which is loop 
suppressed in the standard model (SM). It is expected that the FCNC process is sensitive 
to new physics beyond the SM. On the contrary, the leptonic decays of charm meson, e.g. 
D(s) — * Iv, are tree dominated. New physics in these processes, if exists, should be very 
small. However, recent measurements on the leptonic decays of show that their ratios 
are larger than expectation. In |2|, the authors reviewed the experimental and theoretical 
status of the decay constants of D and Dg. For most approaches to calculate the decay 
constant of D^, theory predictions are smaller than the experiment. In particular, there is a 
3 standard deviation between the lattice QCD calculation which claims a precise prediction 
[sl and the experimental data. A recent updated QCD sum rules analysis provide an upper 
bound which does not reduce the tension between theory and experiment 0]. The above 
discrepancy is sometimes called the '''fos puzzle" . 

To fully understand the fr,^ puzzle, it requires an accurate knowledge of strong interaction 
which is difficult to obtain at present. In this study, we explore the problem from a phe- 
nomenological point of view. Although the treatment has the disadvantage that theoretical 
errors are not well under control, it permits an analytical treatment and its validity can be 
tested by many processes. Our method is a light-front approach jsl, IgI. This is a relativistic 
quark model and its essential element is hadron's light-front wave function. As a successful 
phenomenological model, it has been widely applied to calculate many different meson decay 
constants and form factors. Within this approach, the fn^ puzzl problem really exists. In a 
previous result [sl, the decay constant of Dg is 230 MeV, which is quite different from the 
experimental value of about 270 MeV. How to explain the puzzle of decay constant within 
the light-front approach is our first task. 

Another possible mechanism to explain the decay constant puzzle is to introduce new 
physics effects beyond the SM. There have been many scenarios proposed, e.g. charged 
Higgs and/or leptoquark 0, 0, 0, ll^l- In order to explain the experiment, the new physics 
effects must interfere constructively with the dominant SM contribution. The charged Higgs 
model in [3] is excluded due to its destructive effect. Here, we suggest a new physics scenario, 
unparticle physics [ill, 12| . The unparticle is a scale- invariant stuff with no fixed mass. The 
scale dimension of the unparticle is in general fractional rather than an integral number. 
Many unusual phenomena caused by unparticle are explored 0, Q . The non-trivial phase 
in unparticle theory could induce constructive interference and even a sizable CP violation. 



Some unparticle physics effects in the leptonic decays of B ^ lu are studied in 15|, [iGj. A 
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detailed analysis of D(^s) decays in unparticle physics will be provided in this study. 

If there exits new physics in subprocess c slv transitions, it should also occur in semi- 
leptonic decays of charm mesons, such as -D ^ Klv, Dg — > r],ri' etc.. In most literature, 
the semileptonic decays are either less considered or studied separately. Simultaneously 
studying the charm leptonic and semileptonic decays may help to discriminate different new 
physics models. At the same time, to test the validity of light-front approach, study of the 
semileptonic decays is necessary. 

The paper is organized as follows: In Sec. II, the leptonic and semileptonic decays of D 
and Dg mesons in SM are given. The decay constants and the form factors are calculated 
within the light-front approach. In Section III, we give an analysis of leptonic decays in 
unparticle physics, concentring on the unparticle effects on branching ratios and direct CP 
violation. The last section is devoted to discussions and conclusions. 



II. LEPTONIC AND SEMILEPTONIC DECAYS OF D AND D. MESONS IN SM 



A. L>+x l+u decays in SM 




FIG. 1: The lowest order diagram for the decay of — > I'^vi in SM. 

In SM, the purely leptonic decay of -D(s) Iv occurs via annihilation of quark pair to 
a charged lepton and neutrino through exchange of a virtual W boson. The lowest order 
contribution is a tree diagram which is depicted in Fig. [T]for the decay of l^vi- The 

effective Hamiltonian for subprocess c qlv transition at quark level is 

= ^V,,{qc)v.A{^l)y.A. (1) 

where q denotes d for and s for meson, respectively, and V — A = 7^(1 — 75). 
The weak radiative corrections are so small that they can be safely neglected. The strong 
interactions between c quark and antiquark q by exchange of gluons are incorporated in the 
definition of the decay constant as 

(0|g7M75c|Z^(p)) = tfoP^. (2) 
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Here D represents or to simplify the representation. 
The decay width oi D ^ Iv is obtained straightforwardly by 



T'^{D lu) = ^M^/^M^mf {^~Ml}- 

Because of helicity suppression for spin-0 particle decaying into two spin-1/2 fermions, the 
decay rate is proportional to the lepton mass square mf. 

It seems that the leptonic decays of lu are theoretically simple and clean in 

SM. The weak interaction is well determined and the information of strong interaction is 
encoded in terms of decay constants of meson. If the CKM parameters are known, the 
decay constants of fc^^^ can be extracted from the measured decay ratios by using the above 
equation. From the recent experimental data Br{D^ — > ^Ji^v^) = (4.4 ± 0.7) x 10~^ and 
Br{D+ fi+u^) = (6.2 ± 0.6) x lO'^ (itI, we obtain 

f^P = 221 ± 17 MeV, /^^ = 270 ± 13 MeV. (4) 

where the CKM parameters Vcd = —0.2257 and Vcs = 0.9737 are used. 

Recently, HPQCD and UKQCD collaborations improved their techniques of lattice QCD 
(LQCD). They give a precise prediction for decay constants of D and Ds by 

^LQCD = 207 ± 4 MeV, f^^^^ = 241 ± 3 MeV. (5) 

The most impressive thing about the above results is that the theory error is very small, only 
2% or even smaller for fo^- Comparing Eqs. (jlj) and ([5]), the experiment is consistent with 
theory within the errors for while the experiment data differ from theory prediction by 
about 3a deviations for f^^ ■ Since the s quark in Dg is heavier than the d quark in D meson, 
the calculation for /^^ is expected to be more reliable than fo. The above discrepancy leads 
to so called puzzle. 

At first, we study whether the discrepancy between theory and experiment can be reduced 
within the framework of SM. Our approach is a covariant light-front quark model (LFQM) 
[sl, 0]. This is a relativistic quark model in which a consistent and fully relativistic treatment 
of quark spins and the center-of-mass motion can be carried out. This model has many 
advantages. For example, the light-front wave function is manifestly Lorentz invariant as it 
is expressed in terms of the internal momentum fraction variables which is independent of 



the total hadron momentum. Some applications of this approach can be found in 18|. 
In the LFQM, the decay constant of a pseudoscalar meson is represented by 



fp = ^-ir I dx(fkx^^==(t)p{x,ki_). (6) 
Stt'^ J + k\ 

where A^c = 3 is the color number of QCD, A = uiix + m2(l — x) and mi, 2 represent 
masses of constitute quark and antiquark in the meson. The variable x denotes the light- 
front momentum fraction and k± denotes the intrinsic transverse momentum of the quark. 
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FIG. 2: The decay constants fo and fu^ vs the parameter p. The dashed Hne represents the 
central value and the gray region the error of the experimental data. 



TABLE I: The decay constants of D and Dg in the light-front quark model (in units of MeV). 







Model I 


Model II 


Exp. 


fD 


(MeV) 


200 


221 


221 ± 17 




(MeV) 


230 


270 


270 ± 13 



(f)p{x, k±) is the hadron light-front wave function. In phenomenology, a Gaussian-type wave 
function is widely chosen as 

^(..M^iV^f exp(-^I±S). (7) 

where the normalization constant = 4(-^)^/'^ and ^ = ^^^^^^^^^^ with = mf + k'j_ + 

and Mo = y ^ H — ^ ^ ^ . The parameter (3 in the wave function determines the con- 
finement scale and is expected to be of order Aqcd- The quark masses and /3 are the only 
required phenomenological parameters which makes LFQM predictive. 

In this study, the constitute quark masses are chosen as in rrid = 0.26 GeV and 
rric = 1.40 GeV. So the decay constant depends only on parameter /5. Figur^ plots the 
decay constants of fo and /^^ that vary with p. We can see that the decay constant 
increases as (3 increases, and the relation is nearly a linear function. The slopes df/djS for 
D and Dg are the same to a good accuracy. It is noted that we observe this linear relation 
for the first time although the reason for the relation is unknown. By tuning the parameter 
/?, it is easy to make theory be consistent with experiment. The j3 is within a reasonable 
parameter region and its variation is small, less than 10%. For example, (3ds = 0.59 ± 0.03 
GeV. In fact, as we will show for the semileptonic decays, the recent accurate data provide 
better predictions and determinations of parameters than before. 

Here, we present two different results for decay constants and branching ratios for semilep- 
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TABLE II: The branching ratios for the leptonic decays of D and Ds in the hght-front quark model. 





mode 


Model I 


Model II 


Exp. 


Ds- 




4.5 X 10-3 


6.2 X 10-3 


(6.2 ±0.6) X 10-3 


Ds- 


■> TV 


4.4 X 10"2 


6.0 X 10^2 


(6.6 ±0.6) X 10^2 


Ds - 


-> ev 


1.1 X lO"'^ 


1.5 X lO"'^ 


< 1.3 X 10-"^ 


D - 




3.6 X 10"'' 


4.4 X 10^"^ 


(4.4 ±0.7) X 10-4 


D - 


> TV 


9.6 X 10^"^ 


1.2 X 10^3 


< 2.1 X 10-3 


D - 


> ev 


8.5 X 10"^ 


1.0 X 10"^ 


< 2.4 X 10"^ 



tonic decays in Table HI and HTl The Model I refers to choosing parameters as in [61]. The decay 
constants are chosen as jo = 200 MeV and = 230 MeV. The corresponding parameters 
are fixed to Pd = 0.448 GeV, Pd^ = 0.492 GeV. Obviously, the decay constants in this 
model are smaller than experimental data. In Model II, we make the decay constants fit the 
experiment. In order to fulfill this, the parameters change to Pd = 0.499 GeV, Pd^ = 0.592 



GeV. The experimental data are taken from PDG08 [17|]. The ratio of D tu is predicted 
to be 1.2 X 10^3 in Model II, which is close to the present experimental upper limit. This 
process should be observed soon. For the decays to electron, the ratios are predicted to be 
of order lO^'' or 10^^ which makes them difficult to observe. 

B. Semileptonic decays of D ^ K^*^lv and Ds (piViV') SM 

The semileptonic decays are more complicated in strong dynamics than leptonic processes 
because more hadrons are participating. They play an important role in testing consistency 
of LFQM approach and new physics scenarios. After determining parameters P from the 
leptonic decays, we are able to give predictions for semileptonic decays. In general, the 
variations of P change the charm meson wave function and modify the hadron transitions. 
The chosen processes have the same subprocess c slv transition as the leptonic decays 
Dg — > Iv since we are interested in the "/^i^ puzzle". Thus, the processes to be considered 
include D — > K^*Hu and Dg — > ^(r], r/') decays. They are classified into two categories: 
D — * Plu and D —>■ Vlv depending on whether the final meson is pseudoscalar or vector. 

For semileptonic decay of D meson to a pseudoscalar, i.e. -D(-P) — P(P')/z/, the differen- 



tial partial width is given by [19 1 



rIV n"^ \V P r? 

|_(Z, ^ PI,) = W^^Fr(f )\\ (8) 
where q = P — P' is the momentum transfer and is the invariant mass of the lepton- 
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neutrino pair; p is the final meson momentum in the D rest frame with 



P = \P'\ 



(^Ml - (M - y^y^ (^Ml - (M + 



2M, 



(9) 



D 



where M denotes the final meson mass. Neglecting the light lepton mass, the differential 
partial width is governed by one form factor Fi[q^). The D —>■ P transition form factors are 
defined by 



{P{P')\s^,c\D{P)) = F,{q') 
Then, the total width is 



(p + - — 



+ ^-2 Fo{q')q,. 10 



r{D Piu) 



dV 



dq^' 



For D{P) —>■ K(P')Z+z/ decays where V represents a vector meson, the differential partial 
width is given by 



D 



where p is the vector meson momentum in the D rest frame; the helicity amplitudes H.i{q^) 
are given by the combinations of form factors 



H^{q^) = {MD + Mv)Ar{q^)T 



Mpp 
Md + M^ 



1 



2Mv^q^ 

The D ^ V form factors are defined by 

1 



{Ml -Ml- q^){MD + Mv)Ai{q^ 



2 1\J2^ 



p^M\ 



Md + M^ 



-A2{q') . (13) 



{V{P\e)\-s^,c\D{P)) 



MD + M^ 



V 



{V{P',e)\s^,^,c\D{P)) = z\{MD + Mv)e*Mq 



e* ■ (P + P') 

Md + Mv 



iP + P'),A,iq') 



-2M^ 



e* ■{P + P'] 



V- 



qMsil") - Mq')] 



(14) 



where is the polarization vector of the vector meson V and it satisfies e ■ P' = 0. 

For P ^ P and P ^ V transition form factors, the detailed formulas in the covariant 
light-front approach are given in j^, 0]. We will not display their explicit forms here for 
simplicity. Besides some parameters mentioned in the previous subsection, other necessary 
parameters are: = 0.37 GeV, = 0.3864 GeV, (3k' = 0.2727 GeV, and 0^ = 0.3070 
GeV. They are all taken from 
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For 1] and rj' mesons, the case is complicated by their mixing, i.e., the flavor eigenstates 
are not the physical states. Following |20!], the i] — rj' mixing is given by 

(15) 

where is the mixing angle; r]g = and rfs = ss. From the Feynman diagram, we know 

that only ss contributes to the flnal meson rj^r]') in Dg — > ri^r]'). We need to know and P^, 
in order to calculate the relevant form factors. In the light-front quark model, the parameter 
f3 is extracted from the decay constant. The decay constants are taken to be = —113 





MeV and /*/ = 141 MeV j2l|. The mixing angle is chosen as an averaged value: = 39.3° 
20( 1 . From these above parameters, we obtain = 0.4059 GeV and P^, = 0.4256 GeV. 
In the covariant light-front approach, the formulas of form factors are derived in the frame 
= with = —q'j_ < 0, only values of the form factors in the spacelike momentum region 
can be obtained. The advantage of this choice is that the so-called Z-graph contribution 
arising from the nonvalence quarks vanishes. In order to obtain the physical form factors, an 
analytical extension from the spacelike region to the timelike region is required. The form 
factors in the spacelike region can be parametrized in a three-parameter form as 

i-<'(^)+''(^) 

where F represents the form factors Fi, Ai, and V\ -F(O) represents form factor at 

= 0. The parameters F{0) and a, b are flxed by performing a three-parameter flt to the 
form factors in the spacelike region which can be calculated. We then use these parameters 
to determine the physical form factors in the timelike region. The parameters of a, b and 
F(0) are fltted from the form factors at momentum region —15 GeV^ < g^ < 0. 

The fltted values of F{0) and a, b for different form factors Fi, Ai, A2 and V are given 
in Table UTTl Because Aq and A^ don't appear in Eqs. f|T2l [T3|) . we will not include them 
in Table IIIII Our results in Model I are consistent with those in [6| . The form factors for 
Fi, Ai, A2 in Model I and Model II are nearly equal. There is only a 5-10% difference for 
form factors V(g^) in Models I and II. 

Now, we give predictions for branching ratios of semileptonic decays of D and Ds. The 
results are displayed in Table HVl The experimental data about Dg decays are taken from 
a most recent measurement from CLEO collaborations [22^. About the numerical results, 
some comments are given below: 

(1) The predictions in Model II are in general smaller than those in Model I. For most 
processes, the results in Model II are closer to the experimental data. In other words, the 
semileptonic decays prefer larger decay constants of D and Dg which is indicated by leptonic 
processes. 
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TABLE III: The Form Factors of D and Ds in the hght-front quark model. 





Model I 


Model II 


z? 


F{0} 


a 


b 


F(0) 


a 


b 


^1 


U.79 


1.18 


0.27 


0.78 


1.15 


n o /I 

0.24 


aDK* 


U.DO 


U.o5 


0.03 


0.64 


0.49 


n no 

0.02 


aDK* 
^2 


U.57 


1.U7 


0.33 


0.57 


1 n c 

1.05 


0.27 


V 


U.yo 


1.35 


n /in 

0.49 


0.90 


i.z8 


n /in 

0.40 


T^DsK 


U.7z 


1.27 


0.37 


0.70 


1 1 o 

1.18 


n oo 

0.28 


aDsK* 
^1 


U.OD 


U.d7 


0.09 


0.53 


n c c 

0.55 


n n /I 

0.04 


aDsK* 
^2 




1 1 zL 

1.14: 


u.ou 


u.ou 


1 ns 

i.UO 


U.OO 




0.89 


1.49 


0.76 


0.81 


1.34 


0.51 


pDsV 


0.50 


1.17 


0.34 


0.48 


1.11 


0.25 


pDsV' 


0.62 


1.14 


0.31 


0.60 


1.08 


0.23 


j^Ds<f> 


0.65 


0.60 


0.05 


0.62 


0.49 


0.02 


^2 


0.57 


1.04 


0.37 


0.58 


0.99 


0.26 




1.03 


1.35 


0.57 


0.94 


1.22 


0.39 



TABLE IV: The branching 


ratios of semileptonic 


decays of D and Dg 


in the light-front quark 


model (in units of %). 








decay mode 


Model I 


Model II 


Exp. 




3.90 


3.81 


3.58 ± 0.06 fl7] 


DO ^ i^*-e+z/e 


2.57 


2.38 


2.38 ±0.16 fl7] 




9.96 


9.74 


8.6 ±0.5 [17] 


^ E:*Oe+z/e 


6.50 


6.02 


3.66 ±0.21 [17] 




2.42 


2.25 


2.48 ±0.29 ±0.13 [22] 


D+ ^ r]'e+i^e 


0.95 


0.91 


0.91 ±0.33 ±0.05 [22] 


D+ ^ (/)e+z/e 


2.95 


2.58 


2.29 ±0.37 ±0.11 [22| 



(2) For rj{rj')e^Ve, 4>e^i^e and K*^ e^u^ decays, the theory is in good 
agreement with the experiment. 

(3) It is difficult to understand the decay — > K*^e^Ve where the theory prediction is 
larger than the experiment. Under the isospin symmetry, ^}^^K*-^e+u \ ~ ~ 2.54. But 

the experiment result is BriDO^K*-'e^'u\ ~ ~ 1.54. The discrepancy between the theory 
and experiment may be related to the old puzzle about life time difference between and 
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(4) Our results favor a large rj — t]' mixing angle = 39°. This may be due to our 
neglecting the glue component in rj'. 



III. LEPTONIC AND SEMILEPTONIC DECAYS OF D AND D, MESONS IN 
UNPARTICLE PHYSICS 

A. Leptonic decays in unparticle physics 

As discussed in the introduction, new physics effects must interfere constructively with 
the SM contribution and enhance the rate of leptonic decay. Unparticle physics can provide 
such interference due to the nontrivial phase effects. The scale dimension du of the unparticle 
is in general fractional rather than an integral number. The fractional dimension induce a 
phase factor e"*'^"'^ in the propagator of the unparticle field. 

The scale invariant unparticle fields emerge below an energy scale Ay which is at the 
order of TeV. The unparticle has some peculiar characteristics that make it different from 
the ordinary particle. The interactions between the unparticle and the SM particles are 
described in the framework of low energy effective theory. For our purpose, the coupling of a 
scalar unparticle to two SM fermions (quarks or leptons) is given by an effective interaction 
as 



where On denotes the scalar unparticle fields. The C^f are dimensionless coefficients and 
they depend on different flavors in general. There are two reasons that we don't consider the 
vector unparticle. (1) The transverse vector unparticle does not contribute to the leptonic 
decay of a pseudoscalar meson [15]. (2) Even if there is a nontransverse vector contribution, 
another constraint will suppress it significantly. For the vector unparticle, it is pointed out 
that conformal symmetry puts a lower bound on its scale dimension du > 3 [23^]. If we take 
this constraint seriously, the vector unparticle effects in most processes will be very small 
and are negligible. 

In this study, we are only interested in the effects of the virtual unparticle, thus it only 
appears as a propagator with momentum P and scale dimension dy. The propagator for the 
scalar unparticle field in the timelike momentum region with > is given by [12I, [3] 



■u _ 

•eflf — 




f'l,il-7,)fd''0u + h.c.. 



(17) 




/ 



iP-x 



{Q\T0u{x)0u{m) 



I 





2 sin(rfiYvr) (P^ + ief-du 



where 



Ad, 



167r^/^ T{du + 1/2) 
{2TxYduT{du-l)T{2dur 



(19) 
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The function sinlduir) in the denominator imphes that the scale dimension du cannot be 
integral for > 1 in order to avoid singularity (for djj = 1 the singularity is canceled by 
T{dK — 1) term in A^^). The phase factor e~*^"'^ provides a CP conserving phase which 
produces peculiar interference effects in high energy scattering processes, CP violation in B 
and D decays, etc.. There may be scale symmetry violation after the spontaneous breaking 
due to the coupling of the unparticle to the Higgs [2^. The estimate of the violation is 
difficult, so we neglect it in this study. 




FIG. 3: The lowest order diagram for the decay of Df l^ui in unparticle physics. The double 
dashed lines represent the unparticle. 

The Feynman diagram of unparticle contribution is obtained by replacing the W boson 
to unparticle field. Fig. [3] depicts the lowest order diagram for the decay of l^ui in 

unparticle physics. The amplitude for quark level transition of c ^ qlv is 

where = M|, in D —>■ lu decays and Cq = CcqCiv In the SM, the W boson can 
be integrated out and the interaction of four fermions becomes a local interaction at low 
energy. Because unparticle is different from a heavy particle with a fixed mass, the unparticle 
propagation is a non-local interaction. But is a constant, we can still consider the c qlv 
transition as an effective interaction. Considering Eq. ([6]), the P>^P'^ term in Eq. (120|) can be 
replaced by P'^g^'^ (note that they are not equal) in the final result. The c qlv transition 
is rewritten by 

^^ff = 7^^Tt-^%T i^T e-^'"'^ {<1c)v-a{vI)v-a = re-^^ KT, (21) 
2 sm(dw7r) V / 



where r and are 



A. Cq fMlV" 



2 sin(dwvr) \ A^J CpV^q ' 
kA = dull. (22) 
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1 1.1 1.2 1.3 1.4 1.5 1.6 

du 



FIG. 4: The branching ratio of Dg fii/j^i vs the scale dimension du. The Cg are taken to be three 
values: 1, 10, 41. 

Eq. (pTl) shows a clear physical meaning: the unparticle effects are equivalent to multiplying 
a constant factor with a CP-conserving phase to the SM contribution. 

Combining the SM and unparticle contributions, we obtain the decay rate of D ^ lu 
decays as 

r{D lu) = r^^{D re-''^? . (23) 



Our result is principally consistent with the formulation for B decays in [isl, E^ . 

thors point out a novel CP asymmetry in the leptonic decays caused by 



In 



the CP-even phase of unparticle. The phase (^u mimics the strong interaction phase caused 
by final state interactions. This phenomenon distinguishes unparticle model from other new 
physics scenarios. If there is CP asymmetry m. D ^ Iv decays observed in experiment, it 
would be a clear signal of unparticle physics. Unlike the decay where the CKM parameter 
Vuh contains a CP violating weak phase, the Vcd and V^s in D decays have no weak phase 
in SM. In order to produce CP asymmetry which requires both CP-even and CP-odd phase 
differences, we make an assumption that the coupling coefficient CcqCi^ is complex and 
contains a CP-odd phase (pu, even if its origin is unknown, and then CcqCi^ CcqCi^e~^'^'^ . 
After this assumption, the direct CP asymmetry in Dg —>■ Iv decay is 

V{D~ ^ l-v) - V{Dt - l+v) 



Acp{Dg lu) 



V{D- ^l-u) + V{Dt ^l+u) 
2rsin0i^sin0^ 



(24) 



1 + -f 2rcos(f)ucos(f)w 

The CP violation is caused by the interference between SM and unparticle contributions. 

Now, we discuss the numerical results. We fix the scale parameter An = 1 TeV. The 
coupling coefficients have been defined to be Cg = CcgCiy and Cd = CcdCiu- The SM 
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1 1.1 1.2 1.3 1.4 1.5 1.6 



FIG. 5: The direct CP asymmetry of Dg fiv^ vs the scale dimension cIk. The Cg are taken to 
be: 1, 10, 41. 

parameters are taken from Model I: fo = 200 MeV, fo, = 230 MeV. We don't use Model II 
since new physics is not necessary in it. At first, we give the results for Dg — > ni^ij, where new 
physics effect is expected to be most important. FigJH plots the dependence of branching 
ratio on the scale dimension du at different values of Cg = 1, 10, 41. We gives the results in 
the range 1 < du < 1.6. In order to explain the experiment, the Cg needs to be larger than 
10. For the CP asymmetry, we consider a maximal case, i.e. the CP-odd phase 0^ = 7i/2. 
Fig. [5] plots the dependence of direct CP asymmetry on the scale dimension dy also at 
values of Cg = 1, 10, 41. It is seen that the maximal Acp can reach 35%. Considering 
the experimental constraint for branching ratio, Acp reaches 10% for Cg = 10 and 30% 
for Cg = 44. Thus, the order 10% CP asymmetry is possible in unparticle physics. Our 
predictions for direct CP violation seems large. This is because we adopt a maximal case for 
the CP-odd phase. This phase is unknown and other choices will decrease the predictions. A 
recent measurement from CLEO Collaboration shows no CP violation, Acp = (4.8 ± 6.1)% 



25( 1 . But it does not exclude the unparticle scenario because the experimental errors are 
large and we consider the maximal CP violation in theory. In fact, even a 1% CP violation 
is a support for unparticle theory. Similar results for D /xz/^ decay are given in Figs. [6] 
and El 

The experiments provide stringent constraints on the unparticle parameters. Without 
loss of generality, scale dimension is chosen to be (i^^ = 1.1. Then the observed branching 
ratios of Dg{D) Iv can be used to constrain the coupling coefficients Cd and Cg. Table IVl 
lists the constraints on the Cg and Cd from Dg{D) Iv decays. 



13 





FIG. 7: The direct CP asymmetry of D — fiu^ vs the scale dimension du. The Cd are taken to be: 
0.1, 1, 2.5, 10. 

B. Semileptonic decays in unparticle physics 

For the semileptonic decays, the quark level transition of c — > qlu is nearly the same 
as that in the leptonic decay except the momentum of the unparticle is not equal to that 
of D meson. The unparticle momentum is equal to the lepton pair momentum q, and the 
amplitude of the subprocess is 

By using the equation of motion, v'^^i^ — 75)/ g'^ = in the zero lepton mass limit. Thus, 
a conclusion is obtained: the scalar unparticle contribution is helicity suppressed and van- 
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TABLE V: The constraints on the Cs and Cd from Ds{D) lu decays with du = 1.1. 





Ds fj-v 


Ds TV 


Ds — > ev 


Cs 


41 


46 






D ^ Hu 


D ^TV 


D ^ev 


Cd 


2.5 







ishes for semileptonically decaying to the light lepton (e, /i). In the leptonic decay case, 
because the leading SM contribution suffers the helicity suppression, the unparticle effects, 
although suppressed, play an important role. While for the semileptonic decay, the leading 
SM contribution is not suppressed. So the unpaticle effects are negligible due to helicity 
suppression and the weak coupling with SM particles. The vector unparticle may contribute 
to the semileptonic decays, but it does not enhance the ratio of the leptonic decay and cannot 
solve the fn^ puzzle. 



IV. DISCUSSIONS AND CONCLUSIONS 



We have studied the leptonic and semileptonic decays of D and Ds within SM utilizing a 
light-front quark model. We find that it is not difficult to solve the fo^ puzzle by adjusting 
parameter reasonably. The predictions for semileptonic decays are consistent with exper- 
iment. Although the numerical results depend on the model and the theory uncertainties 
are not under control, the conclusion may be general and model independent. There is a 
sufficient space due to strong interaction uncertainties which can explain the discrepancy 
between theory and experiment. This conclusion is different from claims from lattice QCD 
and QCD sum rules. 

We also study the leptonic decays in unparticle physics. The unparticle induces construc- 
tive interference effects which can enhance the theory to be consistent with the experiment. 
Production of CP asymmetry at percent level is possible. This would be a clear signal for 
unparticle physics. We hope the future experiment can test its validity. 

The discrepancy between the theory and experiment becomes smaller if we use the most 
recent measurement from CLEO collaboration with fu^ = 259.5 ± 6.6 ±3.1 MeV j25|. 
The solution of the fo^ puzzle requires more accurate measurements on the leptonic and 
semileptonic decays of charm mesons. The future BESIII will provide precise determination 
of Ds decay constant with errors to 1% 



26|. 



In conclusion, there is space in SM to interpret the fo^ puzzle without contradictions with 
the other experiments. The observation of CP violation at percent level may be an ideal test 
of the unparticle scenario. 
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